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Abstract 

We use 2-color QCD as a model to study the effects of simultaneous presence of chemical 
potentials for isospin charge, /i/, and for baryon number, \xq. We determine the phase diagrams for 
2 and 4 flavor theories using the method of effective chiral Lagrangians at low densities and weak 
coupling perturbation theory at high densities. We determine the values of various condensates 
and densities as well as the spectrum of excitations as functions of /!/ and fis- A similar analysis 
of QCD with quarks in the adjoint representation is also presented. Our results can be of relevance 
for lattice simulations of these theories. We predict a phase of inhomogeneous condensation (Fulde- 
Ferrel-Larkin-Ovchinnikov phase) in the 2 colour 2 flavor theory, while we do not expect it the 4 
flavor case or in other realizations of QCD with a positive measure. 



1 Introduction 



QCD at non-zero chemical potential has been a subject of numerous studies recently. The phenomenon 
of quark-quark pairing and color-superconductivity (l[ has received renewed attention in view of the 
recent observations [Q] that the superconducting gaps can be large and of possible relevance in astro- 
physics and heavy ion collisions (see (3| for a review). The most interesting finite density phenomena 
(phase transitions, in particular) occur in the non-perturbative regime, thus inviting lattice calcula- 
tion methods. However, the lack of positivity of the Euclidean path integral prevents straightforward 
application of lattice Monte Carlo techniques. QCD at finite density thus remains a theoretical chal- 
lenge. 

This is one of the reasons that certain QCD-like theories with positive Euclidean path integral measure 
have also attracted attention. Examples of such theories are 2-color QCD, QCD with adjoint quarks, 
or 2-flavor QCD with isospin chemical potential ||, [| ^, |8| . In these theories, for each quark with a 
chemical potential \x there is another (conjugate) quark with exactly the same properties, except for 
the opposite sign of /i. The quenched approximation, common in lattice studies, is an approximation 
to such kind of theories, rather than to QCD at finite baryon density with all quarks having equal 
chemical potentials ||. 

In reality, dense baryon or quark matter, such as that arising in the interior of neutron stars or 
heavy ion collisions, is characterized by different, not equal and not opposite, chemical potentials for 
different flavors of quarks. One can describe such situations by the simultaneous presence of chemical 
potentials for baryon charge, [lb, and isospin [ii. The physics of this has been discussed [ 10 1 in the 
context of neutron stars > m) and QCD at large density of isospin and small baryon density 

3> pLs) ||. The most interesting result is the appearance of a phase where ground state breaks 
translational and rotational symmetry. This phase is similar to the so-called FFLO (Fulde-Ferrel- 
Larkin-Ovchinnikov) phase of a BCS superconductor [11]. 

In this paper we investigate the behavior of QCD-like theories, in particular 2-color QCD under 
the influence of both and ^b- We shall use controllable analytical methods in two regimes of these 
theories: a chiral effective Lagrangian, based on the symmetries, in the low density regime, and weak 
coupling perturbation theory in the high density regime. We shall determine the phase diagram in 
the hb-,^1 plane, various condensates, and lowest lying excitations. 

Our goal is to gain an understanding of the physics of the interplay between [ib and The 
theories we study (except N c = Nf = 2) have the advantage of having positive Euclidean path 
integral measure, and can, therefore, be studied on the lattice. Our results can be used as benchmarks 
for numerical lattice studies. 



The presentation of our results is organized as follows. We start with two-colour QCD. In section ^ 
we give the low energy effective Lagrangian upon which the sequel is based. A qualitative argument 
for the structure of the phase diagram is then given in |3| In section |] we analyse the properties of the 
vacuum. As a result we find the classical values for the condensates and densities. Then in section 
[5] we study the expansion about the minimum in order to determine the masses of the low energy 
excitations. The extension of this theory to Nf = 4 with two up and two down flavours is discussed in 
section |6[ The discussion of asymptotically high chemical potentials and the FFLO-phase is in section 
[7|. Finally, in the last section we summarize and comment on the relevance of our results for real QCD. 
The presentation of QCD with quarks in the adjoint representation is placed in the appendix. 
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2 The Effective Theory at Finite \xb and \ij 



Two colour QCD at zero chemical potential is invariant under SU(2iVy) rotations in the chiral limit, 
see e.g [jT^] . This enhanced symmetry (as compared to the SU(/Vy-)xSU(iVj)xU(l) of three colour 
QCD) is manifest in the Lagrangian if we choose to represent this in a basis of quarks ip and conjugate 
quarks rp j|. For Nj = 2 we use 
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where the Pauli matrices T2 and 02 act in colour and spin space respectively. 

The enhanced symmetry manifests itself in the low energy effective theory through the manifold 
of the Goldstone modes associated with the spontaneous breaking of chiral symmetry. In our case 
Nf = 2 and the Goldstone manifold is SU(4)/Sp(4), corresponding to the condensation of \J r \I r - 
SU(4) flavor sextet. The fields on this manifold can be represented by a 4x4 antisymmetric unitary 
matrix E with det E = 1 . 

The effective Lagrangian for the field E of Goldstone modes is determined by the symmetries 
inherited from the microscopic two-colour QCD Lagrangian. In the basis of SU(4) spinors ([l]) the 
mass matrix, baryon charge matrix, and the isospin (third component) charge matrix areQ 
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and the resulting effective Lagrangian is 
F 2 



(2) 



£eff(E) 



-F 2 Ti 



{e{h b B + fnI)^{ti B B + ml) + {ii B B + W /) 2 ) - F 2 mlReTr (ME) . (3) 



The /i-dependent terms in the effective Lagrangian appear through the covariant extension of the 
derivative: 
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d E- 
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{^ B B + niI)E + E(im b B + W /) T 
(fi B B + niI)E + E(n B B + W /) T 



t 



(4) 
(5) 



required by an extended local gauge symmetry [|J|. Therefore, to this order in chiral perturbation 
theory, the Lagrangian at finite [i does not require any extra phenomenological parameters beyond 
the pion decay constant, F, and the chiral condensate in the chiral limit, G. This fact gives predictive 
power to chiral perturbation theory at finite fi. The chiral condensate can be traded for the vacuum 
pion mass using the Gell-Mann-Oaks-Renner relation (m = m u = is the quark mass): 



F 2 m 2 = mG . 



(6) 



1 Note on the notations: We use the normalization of the baryon charge different from ja]. The baryon charge of the 
quark is not 1, as in |j, but 1/2, which comes from l/N c , so that the baryon (diquark in N c = 2) has baryon charge 1. 
For simplicity, we omit subscript 3 from 73. 
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In using the effective Lagrangian constructed above we must of course assume that chiral symmetry 
for N c = Nf = 2 QCD is spontaneously broken. Since we have regarded the hadronic modes as heavy 
the theory is expected to be valid only up to the mass of the lightest non-goldstone hadron. However, 
as we shall show the phase diagram is very rich in the effective range of the theory. 

3 The phase diagram expected from the fij = spectrum 

One can understand main features of the phase diagram in the hb^i plane by considering the low- 
energy spectrum of the theory at finite hb determined in ||. In the vacuum, ^b = fxj = 0, the 
spectrum consists of a degenerate 5-plet: 3 pions ttq, tt^ 1 (B = 0, I = 0, ±1), a baryon and an 
antibaryon q, q* (diquark and antidiquark with B = ±1 and 1 = 0). 

On the horizontal axis, \xi = 0, as a function of fiB, there is a phase transition corresponding to 
the condensation of diquarks. This happens at fiB equal to the diquark mass divided by its baryon 
charge, i.e., at [iB = Tn n . Similarly, at fiB = there should be a transition corresponding to the 
condensation of pions. This happens when fij is equal to the pion mass divided by its isospin charge, 
i.e., at hi = m n . This can also be concluded from the fact that there is a simple discrete symmetry: 
d <-> d accompanied by hb «-> Therefore, the [iB^i phase diagram must be symmetric under 
reflection about fiB = fJ-i line. 

As a function of hb the mass of the pion is constant as long as [iB < rn^, which means the 
transition at hi = m n happens at all such hb- Related to this horizontal line of phase transitions by a 
reflection against the diagonal, there is a vertical line at fiB = m ir, the line where diquarks condense. 
After the two lines meet they must merge into a single \xb = fJ-i line. This is because of the reflection 
symmetry and the fact that the mass of the pion in the diquark condensation phase fis > m n is equal 
to hb according to ||. Two different condensation modes (pion vs diquark) compete at pbB = p>i, thus 
the transition along the \ib = fJ>i line is of first order. 

4 Minimum of the effective Lagrangian 



In order to derive the phase diagram we study the minimum of the classical theory. That is: we seek 
to minimize the static part of the effective Lagrangian 

£ cff (S) = - F 2 Ti {^{hbB + hiI)^ {^bB + nil) + {^bB + nil?) - F 2 m2ReTr (MS) (7) 
of the Lagrangian ([|). Our ansatz for the minimum is 

E = Tim cos a + {Tib cos rj + S/ sin rj) sin a , (8) 
where a and rj are variational parameters, Hm = — A4, while X# and S/ are given by 
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The ansatz is based on the following observations: three terms compete for the alignment of the 
condensate: mAi, hbB, and /U/I. Each of them is independently minimized by Sm, £b, and S/ 
respectively. 
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Inserting the ansatz and using the (anti)commutation relations between the matrices involved we 
reduce the static Lagrangian to 



£ e ff(S) = — F 2 y— (fi B + /jq ) cos 2 q + (fi 2 B — fij) cos(2r/) sin 2 a + fi 2 B + fij + 4m 2 cos aj . (10) 
Extremizing with respect to 77 we find (for 

77 = if fii < fiB 

V = g if M/ > MB • (11) 

Note that for //^ = /i 2 there is degeneracy in the 77-direction. This is a manifestation of the SU(2) 
rotational symmetry (at m = 0) between the flavor which has fi = and its own conjugate field (i.e., 
between d and d in the first quadrant of the fii, /i#-plane). Extremizing with respect to the a-direction 
we find: 

a = if y > 1 
coso; = y if y < 1 (12) 

where 

2 

Mb 

9 

y = —t if fii > ■ (13) 
Mi 

The ansatz (||) is indeed a minimum of the classical theory as we shall prove in the next section. 
Since the vacuum energy is given by the value of the Lagrangian at the minimum we can draw the 
following conclusions 

W) = -%^ = 4Gcosq 



n B = - aC ^ ] = 2F 2 [i B (l + cos (2??)) sin 2 a m = = 2F 2 f i I (l - cos(2t7)) sin 2 a 

aCeff(s) _ / _\ gAgjg) 



= gj l = 4G sin a cos 77 (71") = = 4G sin a sin 77 . 

(14) 

To obtain the diquark and pion condensates we have introduced a diquark source, js, and a pion 
source, ji, in (|3|). The condensates are the derivatives of £ e fj(E) with respect to the sources. The 
introduction of the sources is analogous to that in || and we refer to that paper for details. The 
phase diagram is shown in the left hand side of fig.[l|. It is in agreement with our qualitative argument 
given in section |3|. In particular the densities are discontinuous across at = ^1 > m -n indicating a 
first order phase transition. Note, however, that the vacuum energy, Eq. ([H]), is smooth across the 
Mb = Hi line. The values of the condensates are illustrated in figures Q and 



5 Expanding About the Minimum 

In order to determine the phase diagram we had to find the minimum of the effective Lagrangian 
(||) with respect to S. Having found this, we now determine the masses of low-energy excitations by 
expanding around the minimum E. 
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Figure 1: N c = 2, Nf = 2 (/3d = 1): The figure on the left shows a schematic version of the phase 
diagram. To the right: The ratio of the chiral condensate to its value at a = in the first quadrant 
of the (fJ-B, /i/)-plane. 



5.1 Parameterization of the field S 

A convenient way to parameterize the matrix X in the representation of SU(4)/Sp(4) is based on its 
SU(4) transformation properties [13|: 



(15) 



Here the second equality is imposed by the algebra of the group [13]. The special unitary 4x4 matrix 

U = (16) 
the Lagrangian is minimized by £ = T,m, aligned along 
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With this choice the hermitian field II is restricted by the manifold to [13] 
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(18) 



Where p is real and p and q are complex numbers. The fields p, p, p* are neutral and charged pions 
and q, q* are diquark/antidiquarks. The E will change as a function of fi and we shall discuss the 
parameterization around a generic minimum below. 

The purpose of carrying out an expansion of C e a in the II-fields about the minimum S is twofold. 
First, in order to justify that the ansatz (||) is a minimum, we need to prove that the linear terms in 
II vanish. Second, the expansion of C e g allows us to determine the dispersion relations for the low 
energy excitations. 
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Figure 2: 0y> = 1, Nf = 2: The leftmost plot displays the isospin charge density in the (hb, ///)-plane 
in units of 2F 2 m 7r . On the right hand side is shown the ratio of the pion condensate to the chiral 
condensate at a = 0. The first order phase transition is apparent. 



5.2 Linear Terms 

To expand £ e ff around E defined in (g) we need to choose a parameterization of the IT-field which can 
be used for all values of a and r\. For a = we have already given a parameterization in (jig). For a 
given pair (a > 0, 77) this parameterization of S is no longer valid. We must expand about the rotated 
E according to (^). However, as it is useful not to change the parameterization (|l8|) with a and 77, 
we extract the rotation from U and write it explicitly. To do this we first note that we can write the 
ansatz @ as 

E(a,77) = e -^ S cos,+£™)£ MVjA/ _ = v^EjtfVg,,) . (19) 

Referring to (|l5|) we see that 

S(Q,r?) = U(a,f])^(a,f])U T (a,r]) 

= V {a>v) U(a = 0, »])^, ))) ^ fl )SMVj ))) y t ( afl) [/ T (a = 0, v )V^ v) 

= V {a , v) Z M U 2 (a = 0, V )Vf a>r]) . (20) 

Inserting this in £ e fj @ we find that the form of the Lagrangian does not change provided that we 
substitute rotated values of M. and //b-B + [ill 

^7a,T])^^{oL,r\) = — cos ol — (XI b cos tj + Yjj sin 77) sin ol 
V ( a ,ri)(^BB + fiiI)V^ ar} ) = (hbB cos a - /4b.B£ b Em sin a + m I) cos 77 (21) 

+ {nil cos a — ////E/Ejvf sin a + /ig-B) sin 77 . 

Using this and expanding f/ = exp (OI/2F) to 2 nd order in n we find that the terms liner in II vanish. 
(The cancellation takes place due to the equations (12,13) arising from the constraint (9 a £ st (E) = 0.) 



From this we conclude that the ansatz (||) for the minimum is at least a local extremum of the theory. 
As we now turn to focus on the terms quadratic in II we shall find that this local extremum is indeed 
a minimum. 
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5.3 Dispersion Relations 



The quadratic terms in II of the expansion around S determine the dispersion relations of the low 
energy excitations. For > fij the terms in £ e g(£) at second order in II are 



4ff } (£) = ^Tr {(cUI - ^[1, n]^ - M.B cos a[B, U}5, 
+Ijl 2 b [BT 1 b^'M: n] 2 sin 2 a + m 2 n 2 cos aj . 



(22) 



The corresponding expression for < /i/ is obtained from the above by switching B and / 
symbols. Recall that we have kept the parameterization ( |i~8| ) at the cost of having rotated (^5+ HiI) 
and Ai, see (|2~l"l). Using this parameterization of II we rewrite the above equation as 



4 2 ff } (s) 



(d u p) + (d u p - fj,i5 u0 p)(d u p* + fii5 u0 p*) + (d u q - HB^uoq cos a){d v q* + ^B^uoq* cos a) 



-^f ((« + 1*) 2 + a (pp* + p 2 )) sin2 a + m ^ (p 2 +PP* + ' 



cos a . 



(23) 



This allows us to read off the dispersion relations for each of the 5 modes. 

In the diquark condensation phase a / the q and q* fields mix. The dispersion relations are 
obtained by solving for E = ipo: 
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/x| sin 2 (a)/2 



— (E - fj,B cos a) 2 + p 2 \ 
sin 2 (a)/2 + m 2 cos a 



-(E + \i B cos a) 2 + p 2 
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//^ sin 2 (a)/2 
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(24) 



The mixing between q and q* implies that two of the low energy excitations are linear combinations 
of diquark and anti-diquark modes. We denote them q and q*. As these diquark/anti-diquarks carry 
no isospin it is not surprising that the dispersion relations are the same as found in [p. In the normal 
phase (a = 0) 
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while for a ^ 
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p 2 + /i| ( 1 + 3 cos 2 a) / 2 - {is J^si 1 + 3 cos2 «) 2 / 4 + 4 P 2 cos2 a • ( 26 ) 



Note that one mode is massless for fis > m n . 

There is no mixing in the p, p, p* sector. This implies that we can interpret p, p, and p* as ir°, ir~' 
and tt~ and directly read off the dispersion relations: For a = 
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(27) 
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Figure 3: /3p = 1, Nf = 2: The ratio masses of the pion modes to m n as a function of [iB and [ii- From 
the left m n +, m w o, and m^-. Note that the mass of the tt~ excitation vanishes at /z^ = Hi > n^n- 
The masses of the q* and the q excitations are the mirrors in the hb = ///-plane of m+ and m~ 
respectively. 

while in the diquark condensation phase (a ^ 0) 

7T° E = VP^I 

tt + E = yJp 2 + V 2 B + M 

TT- E = y/pt + nl-M . (28) 

As one might expect we find that the ir + and tt~ modes couple to hi- For /j,b = Hi > rn^ the 7r~ 
becomes massless. For plots of the masses see figure ||. 

The dispersion relations for [ib < Hi 111 the first quadrant are mirrors of the ones above. To obtain 
them one needs only to switch hi ^ hb in (p5|)-(pg|) and charged pions with diquarks/anti-diquarks: 
p^q. _ 

Finally, since the masses of the excitations are positive the value of £ e g(S) with S given in ((^) is 
not only an extremum as shown in the previous section, but is indeed a minimum. We assume that 
this minimum is global. 



6 Positivity and Nf = 4 phase diagram 

Testing these results in lattice Monte Carlo simulations requires care. The fermion determinant, 
though always real in 2-color QCD, is not positive definite if both hb and /// are non-zero. This can 
be seen by factorizing the determinant into a product of the determinants for u and d quarks with 
chemical potentials \i u ^ = Hb i /-*/• Each of these determinants is real but not positive definite if 
Hu,d 7^ 0. So, unless fi u = ±/td, their product is not positive definite. Doubling the number of quarks 
to save positivity will take us to the Nf = 4 theory, with a significantly different phase diagram, which 
we shall now discuss. 

One can understand the phase diagram in this Nf = 4 theory (two up quarks u\ t 2 and two down 
quarks d\ 2) qualitatively using the arguments of sect. |||. The most important difference is that there 

2 We assume that the chiral symmetry is still broken at Nj = 4. The results are trivially extended to 2iV u-quarks 
and 2N d-quarks if chiral symmetry is assumed to be spontaneously broken for that number of flavours. Likewise the 
results on adjoint QCD presented in the appendix may be extended to N u-quarks and TV rf-quarks. 




Figure 4: Phase diagram of the Nf = 4 2-color QCD, determined by a qualitative argument of section 
|6l Solid lines are (second order) phase transitions where certain diquark/antidiquark condensates, 
indicated on the appropriate side of each line, appear. Subscripts refer to flavour. Dashed lines show 
the direction of \x u and \Xd axes. 



are now diquarks which carry isospin: there is an isospin triplet (iii«2, d\d,2, u\d,2 + d\U2) and 3 
isosinglets. As \ib increases, the masses of all 6 diquarks are decreasing as — hb ]!]• Among them 
is u±U2, which becomes the lightest particle with positive isospin. Therefore, for a given fis < m Tr, 
the isospin charge condensation occurs at fij = m n — hb (instead of the horizontal line fij = m n as in 
Nf = 2) - compare figs. [I] and ||. At hb > the diquarks (including the ones carrying isospin) are 
massless and the isospin is condensed already at fii = 0. By the hb ^ Hi reflection symmetry, the 
baryon charge condensation for [ii > m n sets in at fJ,B = 0. The region of the phase diagram to the 
right of (or above) the line Hi = — hb contains a condensate of U\U2, which is favoured by both 
Hb and hi- 

To complete the phase diagram in fig|I| note that the line /i/ = m n — [Xb where u\U2 diquark 
condensation occurs continues also into the negative fii halfplane. Because the theory is invariant 
under [ij — > d <-> u substitution, a line to/ = /xb — m n must exist in the positive \ib^i quadrangle 
(see fig.[|). Below this line the condensate of diquarks dicfo, carrying isospin —1, appears. By [ib ^> Hi, 
d <-> d reflection we obtain another line, \xi = Hb + m n above which the condensate of d\d2 appears. 

An easier way to understand the phase diagram of fig|4] is to consider u and d quarks separately. 
Rotating by 45° and using [i u and [id axes one sees that the diagram is completely symmetric with 
respect to u and d quarks, with condensation of U\U2 and d\d2 diquarks or antidiquarks occurring 
independently once \[i u \ or \fj,d\ exceed the vacuum mass m n . 

We have deferred the rigorous derivation of this phase diagram to future work. However, the phase 
diagram of another theory: 2- (or any-) color QCD with Nf = 2 adjoint quarks (Dyson index /3d = 4) 
is derived in the Appendix. This theory is similar to Nf = 4 (/?d = 1), considered in this section, in 
that it also has diquarks with non-zero isospin: uu and dd. The phase diagrams of these two theories 
are similar. 



10 



7 Large fi and FFLO 



The phase diagram we obtained using chiral perturbation theory is valid only as long as the chemical 
potentials remain small compared to the chiral scale, or m p which is the mass of the lightest hadron 
not included in the chiral Lagrangian. We can, however, study an opposite regime of very large ji, 
when quarks are asymptotically free. 

Let us first consider Nf = 2. On the ub axis, ui = 0, the system is a Fermi liquid of u and 
d quarks with equal chemical potentials. The ground state is a superfluid state with nonzero (ud) 
condensate, which has the same quantum number as the diquark condensate considered previously in 
the framework of the chiral perturbation theory. Thus, it is natural to assume no phase transition on 
the hb axis. Analogously, the ground state at large uj contains ud Cooper pairs, since u and d have the 
same Fermi energy; and there is likely no phase transition separating this phase from the low-/// phase 
of pion condensate. In both cases of large ub arid \ii all fermions acquire a gap A ~ ug~^e~~ c l 9 fl4}| , 
here c = 37r 2 /\/2 and g is the (small) gauge coupling. The proportionality factor has been calculated 
in p]. 

When [ib and [ii are nonvanishing, the chemical potential for the u and d quarks have different 
magnitudes. Generically, {ud) and (ud) condensates cannot be formed unless the mismatch of the 
Fermi momenta is small. Thus, except for small regions near the ub and ui axes, in most of the 
(uB,fJ-i) plane the only condensates that can be formed are (uu) and (dd) (or (dd).) These Cooper 
pairs are color antisymmetric and flavor symmetric, and must hence carry spin or orbital moment and 
break rotational symmetry. 

The region near the diagonal ub = ui requires special consideration. On the diagonal Ud = and 
no d-quarks are present, while the u quarks are paired. Below the scale of the BCS gap for the u 
quark one should expect the u quark to decouple completely. The (E <C A) physics on the diagonal 
is thus Nf = 1 QCD at zero chemical potential. The theory is confining ][l(| and one expects a gap 
for the lightest dd baryon. This baryon cannot be generated unless Ud is larger than this mass gap. 
Therefore, there is a strip along the diagonal where no (dd) or (dd) condensates are present. 

Consider now the regions near the axes, say, the ub one. When ui is small no isospin charge 
is generated. One expects two phase transitions [11, [K|, at ui » 0.71A and 0.75A. In the narrow 
window the ground state is the FFLO phase with spatial varying diquark condensate. Similarly, there 
is a strip of another FFLO phase near the ui axis. We note that the evaluation in [10| applies equally 
well to N c = 3 and N c = 2 and support the occurrence of FFLO in two-colour two-flavour QCD. 

We summarize what is said above in Fig. ||. Since we can solve our model analytically only in 
the two extreme limits, we have no information about how the lines are connected in the region of 
intermediate chemical potentials. Unfortunately, two colour QCD has a sign problem at Nf = 2, so it 
is not clear whether one can find out about it on the lattice. 

For Nf = 4, it is more convenient to work in the variables u u and Ud- When \u u \ ^ \ia<i\, the 
preferable pairing is uu and dd or dd. When Ud = 0, the (E <C A) theory reduces to vacuum QCD 
with Nf = 2 which has confinement and a mass gap for the d\d% diquark. When Ud is less than m^, 
there should not be any d quarks in the system. Thus, the phase diagram should resemble fig. ||, where 
all lines are second order. The u\u<i condensate is always favourable and no FFLO-type transitions 
occur. 

An interesting point is that the non-homogeneous FFLO phase can occur only in the Nf = 2 case, 
and is not expected to arise in the Nf = 4 case. This observation is a special case of a more general 
relation: the FFLO phase in Euclidean QCD is excluded if the measure is positive. This might be 
argued along the following lines, using QCD inequalities [17]. The measure in the QCD-like theories 
is real because of the additional symmetries of the Dirac operator, D. For N c = 2 ((3d = 1) at finite 
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Figure 5: Two colour QCD at small and large chemical potentials: Upper figure shows Nf = 2 and 
lower figure displays Nf = 4. Solid lines are phase transitions and gray areas illustrate regions of 
phase space which remain undetermined. The lines above the clouds are drawn out of scale. Dashed 
lines show the direction of \x u and fid axes. 
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baryon and isospin chemical potential the additional symmetry is 

(T d C^ 5 T 2 Da d C^ 5 T2 = D* with [a d ,a 3 } = and a d = 1 (29) 

where C = ^7072, t~2 acts in colour space, and 03 in flavour space. This relation holds for hb 7^ and 
/i/ 7^ 0. If the measure is positive (by doubling of flavours) the relation leads to a QCD-inequality 
for the diquarks Q. Since o d C^T2 only mix flavours with identical chemical potential so does the 
favored diquark channels for / and [i\ 7^ 0. However, a necessary condition for the FFLO 
phase is that the ud diquark channel is dominant for fig / and [ij 7^ (for degenerate masses). 
Hence positivity excludes the FFLO phase. (The argument for QCD with quarks in the adjoint 
representation, /?d = 4, is equivalent - except in this case we need to assume that quark annihilation 
diagrams can be neglected.) It suggests that a relationship might exist between the presence of a 
translation-invariance-breaking phase and the (absence of) positivity in the corresponding Euclidean 
theory. 

8 Discussion and Conclusions 

One of the main obstacles to the study of QCD at finite baryon chemical potential is the loss of 
positivity of Euclidean path integral measure. Without the positivity standard lattice approach fails. 
One of the ways to gain insight into behavior of QCD at finite density is to study QCD-like theories 
with pseudo-real fermion content. Classified according to the Dyson index of the Dirac operator, /?d, 
such theories are: N c = 2 QCD with fundamental quarks — /?d = 1> and any-7V c QCD with adjoint 
quarks — /?d = 4. 

Such QCD-like theories are studied here and in [||, || ||. We have studied the effects of simultaneous 
baryon and isospin chemical potentials in two limits: low and high density. The low energy limit is 
studied by means of effective low-energy approach - the chiral Lagrangian. Recent first-principle 
lattice studies of two-colour QCD at finite || have confirmed predictions from the effective theory 
at nonzero /ig, with fij = 0. The results of our paper are also testable. For /?d = 1, Nf = 4 and 
f3~E> = 4, Nf = 2 the fermion determinant is positive for all \i& and /x/. 

The low energy effective theory approach, in principle, applies even if positivity of the microscopic 
theory is lost. E.g., it applies to three colour QCD (/?d = 2) at non-zero baryon and isospin chemical 
potential. However, the low-energy effective theory is valid only for energies much less than the 
smallest non-Goldstone hadron mass, and in this regime there is no effect of the baryon chemical 
potential. Low energy effective QCD is, on the contrary, affected if different chemical potentials are 
introduced for different quark flavours |l8|, ^. The isospin case, where two quark chemical potentials 
of the same magnitude, but opposite sign, are introduced, was studied in ||. Comparing that work 
with the one presented here we learn that QCD with (3d = 1, 2, and 4 are very similar on the ///-axis. 
In all cases a pion condensate forms at fii = m n . 

In the HB-, ///-plane the phase diagram distinguishes /3d = 1,2, and 4. For /3 D = 1 the pion 
condensate competes against the diquark condensate, and a novel first order phase transition takes 
place at ^b = l-ii > ni^; for /?d = 2 there is no hb dependence (for hb < ru^)', and for /?d = 4 the fi u 
and n d dependence separate (to lowest order in chiral perturbation theory). 

In the high-density limit we have used weak coupling perturbation theory to study the phase 
diagram. We have found that theories with positivity do not possess a phase of inhomogeneous 
condensation — the FFLO phase. For example, Nf = 2 QCD with adjoint quarks (/?d = 4) is positive 
and does not display FFLO phase. On the contrary, Nf = 2, 2-color QCD (/?d = 1), does have regions 
of FFLO phase at large \ib and /i/ (fiB S> [ii or fif S> hb), while there is no positivity. 
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A Quarks in the adjoint representation 



In this appendix we give results for QCD with quarks in the adjoint colour representation. This 
version of QCD does not have a problem with positivity. For Nf = 1 the determinant is real and posi- 
tive definite (see, e.g., the first reference of §). QCD with quarks in the adjoint colour representation 
belongs to the universality class labelled by Dyson index /?d = 4 jn|, and its low energy limit is also 
described by chiral perturbation theory. Following the line of section ^ to ^ we consider Nf = 2 and 
determine the phase diagram with m = m u = m&. We expect to find a phase diagram similar to that 
of /3d = 1) Nf = 4 since in both cases the lightest excitation on the /i^-axis is a diquark. 



The symmetries of QCD with quarks in the adjoint colour representation allow for a chiral goldstone 
manifold with Nt(2Nf + 1) — 1 degrees of freedom. The goldstone field E may in this case be 
represented by a special unitary symmetric matrix The effective Lagrangian remains to be given 
by (||) provided that we replace M. by 



M = 



( 1 \ 

1 

10 

\ 1 / 



(1) 



The additional degrees of freedom as compared to /3d = 1 {Nf[2Nf + 1) — 1 > Nf(2Nf — 1) — 1)) 
allow for a different vacuum structure. The ansatz for the vacuum field is now 



E 



E^ cos a u + E u sin a u + E^ cos ay + E^ sin a^ 



(2) 



where E^ , E^, E u , and E^ are defined as 



/ 1 \ 
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V o o o o j 



( \ 

1 



y o i o o / 



En 



( \ 
i 

\ i J 

(3) 

With this ansatz the [i u and terms do not mix at the classical level as can be seen from the 
Lagrangian at the minimum 



/ % \ 



i 

\ J 



£ eff (E) = - \-F 2 (- cos(2a u )nl - cos(2a d )/^ + (fi + n\ + 4m^(cos a u + cos a d ) 



(4) 



Extremizing with respect to the a-directions we find: 

,2 



a f = if -£>1 



f = u,d 



cos (a j ) - 



ml if <<x. 



(5) 
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As for /3d = 1 we can prove that £ is a local minimum by expanding around the minimum. The proof 
goes along the lines of the /3d = 1- First we choose a representation of the II field corresponding to 
the direction £ = A4 . Instead of changing this representation as £ rotate with increasing h b and hi 
we again choose to absorb the rotation in hbB + Hi I and Ai. The analogue of Eq. ( |2l| ) is 

V rt (a u ,a d )( / UB-B + ^//)y(a u ,a rf ) = (/j, B B + frl) (£^) 2 cos a u + £ u £^f sin a u 

+ (fi B B + fi I I) (£y) 2 cos ad -l-E^E^ sin a d (6) 
F T (a u ,a d ).A/iy(a n ,a rf ) = cos a n + £ M sin a„ + £}$ cos a d + £ d sin a d . (7) 



with 

V(a u ,a d ) = - iE u e l T^M - iZ u e l ^M . (8) 

Inserting these rotated versions of HbB + hi I and into C e g(E) in Eq. @ we get the /3d = 4 
effective Lagrangian at (a u , a d )- Expanding this to second order in II one finds that the linear terms 
in II drop out due to the extremum conditions of £ e fj (£) in Eq. (||) with respect to a u and ad- Hence 
the ansatz put forward in Eq. (^) is indeed a local extremum. Note that the quadratic terms in II 
mix the u and d sectors, i.e., the u and d sectors only separate at the classical level. As for /3d = 1 we 
will evaluate the condensates and densities classically. Assuming that the local extremum is a global 
minimum we draw the following conclusions 

W>/ = - § %0 1=2Gcosa f > f = u > d 

nB = ~ 9C df^ = 2i?2 [( sin2 a u + sin 2 a d )nB - Hi sin(a u + a d ) siu{a d - a u )\ (9) 

ni = ~ a£ ^j S) = 2F 2 [(sin 2 a u + sin 2 ad) Hi - Hb sm(a u + a d ) sm(a d - a u )] . 

For /3d = 4 and Nf = 2 two of the three diquarks are flavour diagonal, ijr^C^ipf. Introducing diquark 
sources for the diagonal ones we find 

{ipip) f = 2Gsina f , f = u,d. (10) 

The flavour mixing diquark has zero vacuum expectation value outside the hb and Hi axis. On the hb 
axis it is of course degenerate with the two flavour-diagonal diquark condensates. For illustrations see 
figures H and [?]. Finally, let us comment on the possibility of realizing the FFLO state at large hb and 
small 5 hi (or visa versa). The theory allows for two flavour-diagonal diquarks and one flavour-mixing 
diquark. The two flavour diagonal diquark condensates are only mildly affected by 5 hi- The flavour 
mixing state exists only on the hb axis as it competes against the flavour diagonal condensates. The 
FFLO state is therefore not expected to occur. This is consistent with the general relation stated at 
the end of section ^. 
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